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ABSTRACT Viroporins are viral proteins with ion channel (IC) activity that play an 
important role in several processes, including virus replication and pathogenesis. 
While many coronaviruses (CoVs) encode two viroporins, severe acute respiratory 
syndrome CoV (SARS-CoV) encodes three: proteins 3a, E, and 8a. Additionally, pro- 
teins 3a and E have a PDZ-binding motif (PBM), which can potentially bind over 400 
cellular proteins which contain a PDZ domain, making them potentially important 
for the control of cell function. In the present work, a comparative study of the 
functional motifs included within the SARS-CoV viroporins was performed, mostly fo- 
cusing on the roles of the IC and PBM of E and 3a proteins. Our results showed that 
the full-length E and 3a proteins were required for maximal SARS-CoV replication 
and virulence, whereas viroporin 8a had only a minor impact on these activities. A 
virus missing both the E and 3a proteins was not viable, whereas the presence of ei- 
ther protein with a functional PBM restored virus viability. E protein IC activity and 
the presence of its PBM were necessary for virulence in mice. In contrast, the pres- 
ence or absence of the homologous motifs in protein 3a did not influence virus 
pathogenicity. Therefore, dominance of the IC and PBM of protein E over those of 
protein 3a was demonstrated in the induction of pathogenesis in mice. 


IMPORTANCE Collectively, these results demonstrate key roles for the ion channel 
and PBM domains in optimal virus replication and pathogenesis and suggest that 
the viral viroporins and PBMs are suitable targets for antiviral therapy and for muta- 
tion in attenuated SARS-CoV vaccines. 
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= oronaviruses (CoVs) are pathogens responsible for a wide range of existing and 
“ssemerging diseases in humans and domestic or companion animals (1). A CoV 
causing the severe acute respiratory syndrome (SARS-CoV) was identified in Southeast 
China in 2002 and rapidly spread worldwide to more than 30 countries within 
6 months, infecting more than 8,000 people, with mortality in approximately 10% of the 
cases (2, 3). While SARS-CoV has not since reappeared in humans, other CoVs, including 
ones similar to SARS-CoV, are widely disseminated among bats circulating all over the 
world, making future outbreaks possible (4-6). In fact, a novel CoV, the Middle East 
respiratory syndrome coronavirus (MERS-CoV), was identified in September 2012 in two 
human patients with severe respiratory disease in Saudi Arabia (7, 8); since then, the 
WHO has reported 2,144 laboratory-confirmed cases and at least 750 deaths (as of 
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29 March 2018) (http://www.who.int/emergencies/mers-cov/en/). These data indicate 
that emergence of other highly pathogenic CoVs is likely and thus that the study of the 
virus-host interaction is essential to develop antiviral therapies and safe vaccines. 

Viroporins constitute a large class of multifunctional viral proteins with ion channel 
(IC) activity that are widely distributed among different viral families (9); highly patho- 
genic human viruses such as human immunodeficiency virus 1 (HIV-1), hepatitis C virus 
(HCV), influenza A virus (IAV), rotavirus (RV), enterovirus, and CoVs such as SARS-CoV 
and MERS-CoV encode them (10-16). Viroporins promote several steps of the virus 
replication cycle, including entry, genome replication, morphogenesis, and release from 
the infected cell (17, 18). Several viroporins have important roles in viral pathogenesis, 
promoting ion imbalances within cells (13, 19-21) or disrupting cellular pathways 
through protein-protein interactions (22). Given their potential as antiviral targets, there 
is substantial interest in the study of these proteins (18, 23). 

Several CoVs, such as MERS-CoV, HCoV-229E, HCoV-OC43, and porcine epidemic 
diarrhea virus (PEDV), encode two viroporins (24-26), but, remarkably, SARS-CoV en- 
codes three: proteins 3a, E, and 8a (14, 27, 28). The 3a protein is 274 amino acids (aa) 
in length with three transmembrane domains (TMDs). It is the largest SARS-CoV 
accessory protein and is likely involved in virus release (27) and pathogenesis (29). It 
causes membrane rearrangements in infected cells, leading to an increase in levels of 
intracellular vesicles that may facilitate nonlytic release of viral particles (30). Further- 
more, it colocalizes with M protein, which, together with E protein, is essential for virus 
assembly, supporting the notion that the 3a protein is important for SARS-CoV assem- 
bly or budding (31, 32). However, most studies have been based on 3a overexpression 
and little is known about the relevance of this protein in the context of natural 
infection. 

SARS-CoV E protein is an integral membrane protein of 76 aa with only one TMD. 
The residues responsible for E protein IC activity have been previously identified (33, 
34). E protein IC activity is important for SARS-CoV fitness and pathogenesis, since both 
were diminished in its absence (34). 

The 8a protein found in SARS-CoV-infected human cells resulted from a 29- 
nucleotide (nt) deletion in open reading frame 8 (ORF8) that occurred after the virus 
crossed species to infect humans (35). ORF8 genes encode two proteins, ORF8a and 
ORF8b, which represent proteins of 39 and 84 aa, respectively. Overexpression assays 
showed that ORF8a induces apoptosis through a mitochondrion-dependent pathway 
(36) and has IC activity (28). 

Viroporins and cellular IC proteins often rely on protein-protein interactions for 
clustering of ICs at proper locations in the cell (37-39). These interactions are mediated 
between PDZ domains and PDZ-binding motifs (PBMs), peptide sequences that are 
most frequently located at the C terminus of the IC proteins (40, 41). PDZ domains are 
protein recognition sequences, 80 to 90 aa in length, and constitute a large family of 
globular domains found in prokaryotes and eukaryotes. There are more than 400 
cellular protein isoforms containing a PDZ domain in the human proteome (42). The 
PBM core sequence includes 4-aa residues, numbered from the C terminus (p0), which 
is always hydrophobic, to the N terminus (p-1, p-2, and p-3). There are three classes of 
PBMs, depending on the identity of residue p-2: class | for Thr/Ser, class Il for any 
hydrophobic residue, and class Ill for Glu/Asp. Protein-protein interactions involving 
PDZ domains modulate cellular pathways important for viral replication, dissemination 
in the host, and pathogenesis (43). Furthermore, some PDZs also bind PBMs located in 
the internal region of proteins or lipids (44, 45). Of the three SARS-CoV viroporins, both 
proteins 3a and E have a class Il PBM at their C terminus; while the PBM of E protein is 
involved in pathogenesis (46, 47), the role of the PBM of protein 3a, and of similar 
motifs present in other CoV proteins such as MERS-CoV proteins E and 5, remains 
unknown. 

In the present work, using mutational analysis, we showed that only the 3a and E 
proteins were clearly involved in SARS-CoV replication and virulence. Neither single 
deletion of the IC activity or of PBM from protein 3a diminished SARS-CoV replication 
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and virulence, in contrast with E protein, which required both for virulence (34, 46). 
Protein 3a IC activity was characterized in planar lipid bilayers, showing that it forms 
non-voltage-gated ion channels. Furthermore, we identified residues located in TMD2 
and TMD3 of the 3a protein that were involved in its IC activity. The potential 
interdependence of the three viroporins described in SARS-CoV was studied by dele- 
tion of single viroporins or of different combinations of two viroporins. The variant 
missing both the 3a and E proteins was not viable, indicating that the presence of at 
least one of the proteins is essential for virus viability. Furthermore, it was shown that 
either protein should maintain its PBM to compensate for the absence of the other 
full-length protein. These results suggest that PBMs interact with cellular proteins with 
PDZ domains to change cell metabolism, enhancing virus replication or pathogenicity. 
Identification and inhibition of specific cellular pathways affected by these interactions 
may be crucial for the identification of new antiviral strategies. 


RESULTS 

SARS-CoV viroporins E and 3a were required for efficient replication in vitro 
and in vivo. To study the role of the SARS-CoV viroporins in virus replication and 
virulence, three mutant viruses, each lacking one gene (recombinant SARS [rSARS]- 
CoV-MA15-A3a, -AE, and -A8a), were engineered from a mouse-adapted infectious 
cDNA clone (MA15) (48, 49). Analysis of the growth kinetics of each mutant in Vero E6 
cell supernatants (Fig. 1A) was used to determine their requirement for replication. 
Cell-associated virus was also analyzed at 24 and 48 h postinfection (hpi), showing 
results similar to those observed for the released virus (see Fig. S1 in the supplemental 
material). The A3a and AE mutants grew to lower titers than the parental wild-type (wt) 
virus (Fig. 1A). However, while the AE mutant showed 100-fold-lower titers (around 8 x 
10° PFU/ml), A3a titers decreased slightly (3-fold) (3 X 107 PFU/ml). These results show 
that both proteins were required for optimal virus replication in cell culture. In contrast, 
the A8a virus reached peak titers (9 X 107 PFU/ml) similar to those observed for the 
parental virus. 

To evaluate the requirement for protein 3a, E, or 8a for optimal virus growth in vivo, 
BALB/c mice were infected either with rSARS-CoV-MA15 or with each of the viroporin 
deletion mutants SARS-CoV-MA15-A3a, -AE, and -A8a, and viral titers in lungs were 
determined at 2 and 4 days postinfection (dpi) (Fig. 1B). The highest titers were reached 
at 2 dpi, and the titers decreased in all cases by between 5-fold and 40-fold at 4 dpi, 
with the parental and A8a viruses achieving the highest titers in lung tissue (around 
108 PFU/g at 2 dpi and 2 X 107 PFU/g at 4 dpi). However, compared to the parental 
virus, titers were reduced by 1 and 2 log units in the case of A3a and AE virus, 
respectively. Interestingly, the decrease in virus titers after SARS-CoV-MA15-A3a infec- 
tion was greater in vivo than in vitro. Thus, proteins E and 3a were shown to be critical 
for both in vitro and in vivo virus replication. 

SARS-CoV viroporins E and 3a were both associated with virulence in a mouse 
model. To evaluate the relevance of SARS-CoV E, 3a, and 8a viroporins for virulence, 
BALB/c mice were either subjected to mock infection or infected with parental rSARS- 
CoV-MA15 or with one of the deletion mutants rSARS-CoV-MA15-A3a, -AE, and -A8a. 
Clinical disease and survival were monitored through 10 dpi (Fig. 2). Mice infected with 
viruses lacking either E protein or 3a protein recovered from infection with 100% 
survival, although mice infected with the A3a virus showed mild disease symptoms. In 
contrast, mice infected with the parental virus or the A8a virus developed manifesta- 
tions of serious disease (lethargy and ruffled fur) starting from 2 dpi. These mice all died 
by 6 dpi, clearly showing that both E and 3a proteins were involved in SARS-CoV 
virulence in the mouse model, while 8a did not seem to play a major role. 

Characterization of the IC activity of protein 3a in planar lipid bilayers. The IC 
activity of protein E is required for SARS-CoV replication and virulence (34). However, as 
the relevance of the IC activity of the 3a protein was not known, we studied 3a protein 
in planar lipid bilayers and identified the amino acids involved in ion conductance. This 
system was used because its high sensitivity allows the detection of electric currents of 
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FIG 1 Growth kinetics of SARS-CoV viroporin-defective mutants. (A) Subconfluent monolayers of Vero E6 
cells were infected with wild-type (WT) (black filled circles), AE (blue filled squares), A3a (red filled circles), 
and A8a (green filled diamonds) SARS-CoV at a MOI of 0.001. Culture supernatants were collected at 4, 
24, 48, and 72 hpi and titrated by plaque assay. The results are representative of three replicate 
experiments. (B) Groups of six 16-week-old BALB/c mice were infected with 100,000 PFU of either the 
parental virus (WT, black columns) or genetically engineered viruses lacking E protein (AE, blue columns), 
3a protein (A3a, red columns), or 8a protein (A8a, green columns). At 2 and 4 dpi, 3 mice from each group 
were sacrificed to determine lung virus titers. Data summarize two replicate experiments. Data represent 
means + standard deviations (SD). *, P value <0.1; **, P value <0.01; ***, P value <0.001. 


a single ion channel (50). To this end, a baculovirus was engineered to express the 
parental 3a protein in Sf-9 cells. Conductance of purified protein 3a was evaluated in 
the presence of KCI in planar lipid bilayers with a biologically relevant mix of 
1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/1,2-dioleoyl-sn-glycero-3-phospho- 
L-serine (DOPS)/1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE) with ratios of 
3:1:1 (wt/wt), which is a composition similar to that of intracellular organelle mem- 
branes, such as the endoplasmic reticulum (ER)-Golgi intermediate compartment (ER- 
GIC). Single-channel conductance was estimated from a statistical analysis of the 
current jump amplitudes. This procedure allows a reliable estimate of the most prob- 
able value of current change every time a new channel is inserted or disappears. 
Although several channels were being inserted, the magnitude of the current through 
a single channel could be discriminated. Current jumps corresponding to 201 inde- 
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FIG 2 Virulence of SARS-CoV viroporin-defective mutants. Groups of five 16-week-old BALB/c mice were 
subjected to mock infection (PBS, gray filled circles) or infected with 100,000 PFU of either the parental virus (wt, 
black filled circles) or genetically engineered viruses missing E protein (AE, blue filled squares), 3a protein (A3a, 
red filled circles), or 8a protein (A8a, green filled diamonds). Mean levels of weight loss (left graph) and survival 
(right graph) through 10 dpi are shown for each group. Data summarize two replicate experiments with 
equivalent results. Error bars represent the standard deviations of mouse weight data. 


pendent events were measured under conditions of an applied voltage of +100 mV. 
Histograms of the current jump amplitudes of the recorded traces showed that the 
most frequent events corresponded to single-channel conductance of 16 pA (Fig. 3A 
and B). Simultaneous bursts of two or three 3a ion channels were also observed, 
although with much lower frequency (Fig. 3B). 

To test if 3a behaves as a voltage-gated IC, its activity was also measured in planar 
lipid bilayers under conditions of different voltages in the presence of monovalent 
(NaCl) and divalent (CaCl) cations (Fig. 3C). In all cases, a linear current-voltage (l-V) 
relationship was obtained, demonstrating that the channel displayed resistance (ohmic) 
behavior for both positive and negative potential. These results indicate that the 
protein 3a IC was neither open nor closed at specific electric potentials. 

Measurement of the reversal potential (F,.,) of an ion channel, which is defined as 
the voltage that needs to be applied to yield zero electric current when there is an ion 
concentration gradient across the membrane, is the method of choice to quantify ion 
selectivity. Determination of the sign of the Ee, provides a quick estimation of the 
channel selectivity, that is, of its preference for cations or anions (51, 52). By comparing 
the measured E,.,, to the theoretical E e, that would be obtained in the case of a neutral 
pore (i.e., representing only the difference between cation and anion intrinsic mobili- 
ties), the selectivity of the ion channel can be inferred. The theoretical E,.., is calculated 
using the Goldman-Hodgkin-Katz (GHK) equation, replacing the permeability ratio 
P_,/P_ by the solution diffusion coefficient ratio D,/D_ (53). Higher, lower, or similar 
measured E£,,,, values indicate anion selectivity, cation selectivity, or no selectivity, 
respectively. Interestingly, in the presence of monovalent ions (Na+ and K+), the 
protein 3a IC showed weak cation selectivity. However, in the presence of Cat*, the 
channel behaved as a neutral channel with no preference for anions or cations 
(Table 1). Taken together, these results indicated that at least Nat, K+, and Ca++ were 
conducted through the 3a protein IC. 

Identification of amino acids involved in protein 3a IC activity. In order to 
identify the amino acids necessary for protein 3a IC activity, a set of recombinant 
baculoviruses (rBV) expressing mutated 3a proteins was engineered. Amino acid sub- 
stitutions were created to disrupt the IC activity of 3a with minimal impact on its 
three-dimensional structure, mutating residues predicted to face the lumen of the pore. 
As the 3a protein structure has not yet been experimentally determined, in silico 
models were used to select the residues potentially facing the lumen of the pore (54, 
55). As the 3a protein has three TMDs, mutants with changes to a single TMD (TMD1—, 
TMD2~—, or TMD3~) or to two TMDs (TMD[2,3]~) were engineered (Table 2). Mutant 3a 
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FIG 3 Characterization of the SARS-CoV 3a protein ion channel. (A) Recording of a single-channel insertion of SARS-CoV 
3a protein. (B) Histogram of current jump amplitude (right) at +100 mV in 500 mM KCI, composed of values from 201 
recording events. (C) SARS-CoV 3a protein voltage-independent ion channel. The 3a protein showed a linear current-voltage 
relationship. Displayed data correspond to representative l-V plots from reversal potential experiments performed with 
500/50 mM solutions of monovalent (NaCl, red filled squares) and divalent (CaCl, blue filled triangles) cations. Each 
experiment was performed at least three times; the lines represent linear regression fits of data points. 


proteins were expressed in insect cells and purified, and their IC activity was evaluated 
in planar lipid bilayers (Fig. 4). The TMD1~— mutant retained IC activity, but the TMD2-, 
TMD3~—, and TMD[2,3]~ variants did not, consistent with the importance of TMD2 and 
TMD3, as predicted by in silico models (54, 55). 


TABLE 1 Results of protein 3a reversal potential experiments performed with 500/50 mM 
salt solutions 





Ey (mV)2 
lon solution Experimental Reference 
NaCl —19.1 + 14.8 +8.86 
KCI —13.0 + 4.0 +0.73 
CaCl, +18.8 + 6.0 +20.3 





-Experimental reversal potential (F,.,) values represent the averages of results from at least 7 independent 
experiments. Reference Eey values represent theoretical values for a neutral pore. 
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TABLE 2 3a protein ion channel mutations 





Mutant Mutations 

TMD1— S40A, S58A 

TMD2-— Y91A, H93A 

TMD3— Y109A, Y113A, Q116A 

TMD[2,3]~ Y91A, H93A, Y109A, Y113A, Q116A 


To resolve the exact residues necessary for protein 3a IC activity, a complementary 
set of baculoviruses incorporating single amino acid substitutions within mutants 
TMD2- (Y91A and H93A) and TMD3- (Y109A, Y113A, and Q116A) was generated. These 
mutant proteins were expressed and purified, and their IC activity was evaluated. TMD2 
point mutants Y91A and H93A and TMD3 point mutant Y109A completely abrogated 
protein 3a IC activity, whereas TMD3 point mutants Y113 and Q116 showed conduc- 
tance that was equivalent to and only moderately decreased from that seen with the 
wt protein, respectively. Therefore, these results identified 3 amino acids that could be 
mutated to eliminate protein 3a IC activity. 

Protein 3a IC activity was not required for SARS-CoV replication and virulence. 
To study the relevance of protein 3a IC activity in virus replication and virulence, the 
following collection of full-length rSARS-CoVs with and without protein 3a IC activity 
was generated by introducing specific mutations into the 3a gene: rSARS-CoV-MA15- 
3a-TMD1—, -TMD2~, -TMD3~, -TMD[2,3]—, -Y91A, -H93A, -Y109A, -Y113A, and -Q116A. 
All these viruses were similar with respect to growth kinetics in Vero E6 cells (Fig. 5A), 
indicating that replication was not significantly affected by altered protein 3a IC activity. 

The requirement for protein 3a IC activity in vivo was also studied by measuring the 
titers of 3a IC mutants in the lungs of infected BALB/c mice at 2 and 4 dpi (Fig. 5B). Peak 
titers were reached at 2 dpi and had decreased by around 1 log unit at 4 dpi in all cases. 
Every mutant showed replication levels similar to those seen with the wt strain (around 
1 X 107 PFU/g of lung tissue at 2 dpi and 4 X 10° PFU/g at 4 dpi), with the exception 
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FIG 4 Effect of mutations on the ion channel activity of SARS-CoV 3a protein. Recombinant 3a protein 
variants were reconstituted in artificial lipid bilayers, and their IC activity was tested in 500 mM KCI 
solutions. Mean conductance values were measured for variants showing IC activity. Negative control 
(C7) data indicate conductance values obtained in the absence of any protein; error bars represent the 
standard deviations of data obtained in at least 100 independent measurements. 
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FIG 5 Growth kinetics of SARS-CoV mutants targeting 3a protein ion channel activity. (A) Subconfluent 
monolayers of Vero E6 cells were infected at a MOI of 0.001 with wild-type SARS-CoV (WT, black filled 
circles) or with variants with mutations affecting 3a protein TMD1 (TMD1-, light green filled circles), TMD2 
(TMD2-, light blue filled triangles), TMD3 (TMD3-, ochre filled triangles), or both TMD2 and TMD3 
(TMD[2,3]-, light brown filled triangles) or residue Y91 (Y91A, purple filled diamonds), residue H93 (H93A, 
deep blue filled diamonds), residue Y109 (Y109A, orange filled diamonds), residue Y113 (Y113A, dark 
brown filled diamonds), or residue Q116 (Q116A, red filled diamonds). Culture supernatants collected at 
4, 24, 48, and 72 hpi were titrated by plaque assay. Results are representative of three replicate 
experiments. For the sake of clarity, SD data are not shown but the values were, in all cases, lower than 
5%. (B) Groups of six 16-week-old BALB/c mice were infected with 100,000 PFU of either the parental 
virus (WT, black columns), or mutants TMD1~ (light green columns), TMD2~ (light blue columns), TMD3— 
(ochre columns), TMD[2,3]~ (light brown columns), Y91A (purple columns), H93A (deep blue columns), 
Y109A (orange columns), Y113A (dark brown columns), and Q116A (red columns). At 2 and 4 dpi, 3 mice 
from each group were sacrificed to determine virus titers. Data summarize two replicate experiments. 
Data represent means + SD. *, P value <0.1; **, P value <0.01; ***, P value <0.001. 


of TMD3~, which had titers at least 1 log unit lower than those seen with the rest of the 
viruses (3 X 10° PFU/g at 2 dpi and 1 X 10* PFU/g at 4 dpi) (Fig. 5B). These results 
indicated that protein 3a IC activity was not essential for SARS-CoV replication in mouse 
lungs. 

The requirement of protein 3a IC activity for SARS-CoV virulence was studied in two 
independent experiments. In the first, BALB/c mice were infected with rSARS-CoV-MA15 
(virulent virus control), rsARS-CoV-MA15-A3a (attenuated virus control), or one of the 
rSARS-CoV mutants TMD1~, TMD2~, TMD3~, and TMD[2,3]~ (Fig. 6A), with clinical 
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FIG 6 Virulence of SARS-CoV 3a ion channel mutants. Groups of five 16-week-old BALB/c mice were 
subjected to mock infection (PBS, gray filled squares) or infected with 100,000 PFU of the parental virus (WT, 
black hexagons) or with genetically engineered mutants lacking the 3a protein (A3a, dark green filled 
circles), and (A) mutants with altered TMD1 (TMD1-, light green filled circles), TMD2 (TMD2~, light blue 
filled diamonds), TMD3 (TMD3~, ochre filled diamonds), or both TMD2 and TMD3 (TMD[2,3]~, brown filled 
diamonds). (B) Additional comparisons with 3a mutants Y91A (purple filled diamonds), H93A (deep blue 
filled triangles), Y109A (orange filled triangles), Y113A (brown filled triangles), or Q116A (red filled 
diamonds) were performed. All mice were evaluated for weight loss (left) and survival (right) through 10 
dpi. Data summarize two replicate experiments with equivalent results. Error bars represent the standard 
deviations for mouse weight. 


disease and survival evaluated for 10 days. All mice infected with the parental virus or 
the TMD1~ mutant showed disease symptoms starting at 2 dpi, and all died at 5 or 
7 dpi, respectively. Mice infected with the TMD2~ mutant showed acute disease 
starting at 2 dpi, and 80% of the mice had died by between 4 and 7 dpi. In contrast, 
mice infected with either TMD3~ or TMD[2,3]~ mutants recovered from the disease 
with 100% survival, similarly to mice infected with the attenuated SARS-CoV-A3a 
variant (Fig. 6A). During the experiment, viruses were recovered from the lungs of 
moribund mice, and the 3a gene was sequenced. No compensatory mutations restor- 
ing IC activity were identified in any case. The results suggest that the TMD3~ and 
TMD[2,3]~ mutants were attenuated in an IC-independent manner, since the TMD2— 
mutant was only marginally attenuated. 

In the second virulence experiment, mice were infected with rSARS-CoV incorpo- 
rating protein 3a point mutations (Y91A, H93A, Y109A, Y113A, and Q116A), with the 
parental and A3a variants serving as virulent and attenuated controls, respectively. All 
of the 3a point mutants caused severe disease with 100% mortality by 6 dpi, similarly 
to mice infected with the wt parent and in contrast to those infected with the A3a 
variant, which survived (Fig. 6B). No compensatory mutations restoring IC activity were 
identified in this experiment, further confirming that protein 3a IC activity was not 
essential for SARS-CoV virulence in the mouse model used. 
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FIG 7 Requirement of the PBM of SARS-CoV 3a protein for replication and virulence. (A) (Left panel) 
Subconfluent monolayers of Vero E6 cells were infected with wild-type (WT; black filled circles), A3a 
(green filled circles), or 3a-PBM~ (red filled diamonds) SARS-CoV at a MOI of 0.001. Culture supernatants 
collected at 4, 24, 48, and 72 hpi were titrated by plaque assay. Results are representative of three 
replicate experiments. (Right panel) Groups of six 16-week-old BALB/c mice were infected with 
100,000 PFU of the parental virus (WT, black columns) or of a SARS-CoV variant lacking the protein 3a 
PBM (3a-PBM~, red columns). At 2 and 4 dpi, 3 mice from each group were sacrificed to determine virus 
titers. Data summarize two replicate experiments. Data represent means + SD. (B) Groups of five 
16-week-old BALB/c mice were subjected to mock infection (PBS, gray filled squares) or infected with 
100,000 PFU of the parental (wild-type) virus (WT; black filled hexagons) or of genetically engineered 
mutants lacking the 3a protein (A3a, green filled circles) or lacking the protein 3a PBM (3a-PBM~, red 
filled diamonds). Mean levels of weight loss (left graph) and survival (right graph) through 10 dpi are 
represented for each group. Data summarize two replicate experiments with equivalent results. Error bars 
represent the standard deviations of mouse weight data. 


The PBM of protein 3a was not required for SARS-CoV replication and viru- 
lence. To analyze the requirement for the PBM of protein 3a for replication and 
virulence, a virus lacking a functional PBM in the 3a protein (3aPBM_) was engineered. 
As ORF3a partially overlaps ORF3b, the protein 3a PBM core sequence (SVPL) was 
disrupted with amino acid substitutions (GMSM), with codons carefully selected to 
ensure that protein ORF3b was not mutated. Growth of the 3aPBM~ mutant in Vero E6 
cells and in the lungs of infected mice was the same as that seen with rSARS-CoV-MA15 
(Fig. 7A). Also, the 3aPBM~ mutant was as pathogenic as the parental virus (Fig. 7B). 
These results indicated that, in the mouse model, replication and virus virulence were 
independent of the PBM of protein 3a. 

Simultaneous requirement of viroporins by SARS-CoV. In order to study the 
interdependence of SARS-CoV viroporins, all possible combinations of single-, double-, 
and triple-deletion mutants were engineered (Table 3). All combinations were effi- 
ciently rescued, with the remarkable exceptions of the triple mutant [rSARS-CoV-MA15- 
A(3a,E,8a)] and the one lacking both E and 3a proteins [rSARS-CoV-MA15-A(3a,E)]. The 
AE virus and double-deletion mutants A(3a,8a) and A(E,8a) showed significantly re- 
duced titers (1 X 10°, 2 X 10°, and 6 X 10° PFU/ml, respectively) compared to the A3a 
and A8a mutants and the parental virus. Although a role for protein 8a was not 


May/June 2018 Volume 9 Issue 3 e02325-17 


mBio’ 


mbio.asm.org 10 


Hio'wseoiqu Aq pausiqnd - 810Z ‘pz Ae uo Hio"wseolqu wo papeojuMoq 


SARS-CoV Viroporins in Replication and Pathogenesis 


TABLE 3 Simultaneous requirements of viroporins by SARS-CoV 
SARS-CoV viroporin 





Virus name 3a E 8a Viral titer (PFU/ml) 
SARS-CoV wt + + + (4.0 + 1.2) X 107 
SARS-CoV A8a + + - (5.0 + 2.1) X 107 
SARS-CoV A3a - + + (1.0 + 1.9) X 107 
SARS-CoV AE + - + (1.0 + 0.8) x 10° 
SARS-CoV AIE, 8a] + - - (6.6 + 1.4) X 105 
SARS-CoV A[3a, 8a] - + - (2.4 + 1.1) X 106 
SARS-CoV A[3a, E] - - + <2.0 X 1014 
SARS-CoV A[3a, E, 8a] - - - <2.0 X 1014 


“Data were below the detection threshold. 


described for virus replication, the titers of viruses missing 8a and 3a or E were 
significantly lower than the titers of A3a and AE viruses. Also, every one of the viable 
mutant viruses which lacked either 3a or E proteins showed plaque sizes that were 
significantly smaller than those seen with the wt and A8a viruses. 

The results showed that at least either the 3a protein or the E protein must be 
present for virus viability, which indicates that just one of these proteins could provide 
the activities required for virus growth. 

SARS-CoV 3a protein subcellular localization. To begin to understand the nature 
of the replacement of E protein by 3a, or vice versa, in the replication of SARS-CoV, we 
assessed the cellular localization of each protein. The SARS-CoV E protein localizes to 
the ERGIC (56), but less is known about protein 3a localization (27). The potential 
colocalization of the 3a and E proteins was first studied by infecting Vero E6 cells with 
wt rSARS-CoV. Confocal microscopy analysis of rsARS-CoV-MA15 at 24 hpi showed that 
the proteins were located in different cell subcompartments during infection, as the 
Pearson's coefficient value was below 0.6 (see Table S1 in the supplemental material) 
(57); protein 3a mainly accumulated at an unidentified perinuclear compartment 
different from the ERGIC (Fig. 8A). E protein IC activity promotes virulence by releasing 
calcium from the ERGIC, leading to inflammasome activation (34, 58). We showed that 
protein 3a IC activity is not involved in virus virulence. To extend these results, we next 
determined whether 3a was localized to any of the known cellular calcium reservoirs, 
i.e. the ER, the Golgi apparatus, and the mitochondria (59), using specific markers, 
including protein disulfide isomerase (PDI) (ER marker), 58 K (Golgi marker), and 
aconitase 2 (mitochondrion marker) (Fig. 8B). Also, other subcellular compartments 
were analyzed for colocalization with the 3a protein by the use of antibodies (Abs) 
against Na*/K* ATPase (plasma membrane marker), Rab5 (early endosome marker), 
Rab7 (late endosome marker), and LAMP-1 (lysosome marker) (Fig. 9). Pearson’s coef- 
ficient was below 0.6 in all cases (Table S1). Overall, 3a protein did not localize at any 
of the main intracellular calcium storage locations, suggesting that it is not located at 
a site that facilitates increased cytosolic calcium levels. Collectively, these results 
suggest that the IC activity of proteins 3a and E, while shared by the two proteins, is not 
the function responsible for the replacement for the full-length E and 3a proteins. 

Identification of candidate motifs responsible for the replacement between E 
and 3a proteins. To identify which E protein domain was essential in the absence of 
full-length 3a protein, a set of A3a mutants was engineered with substitutions or 
deletions throughout the E protein sequence (Fig. 10A), including (i) alanine substitu- 
tions in the N-terminal domain [rSARS-CoV-MA15-(A3a,EA1)); (ii) disruption of the TMD, 
abrogating E protein IC activity [(A3a,E-N15A)] (34); (iii) short in-frame deletions 
throughout the C-terminal domain [(A3a,EA2), (A3a,EA3), (A3a,EA4), (A3a,EA5), and 
(A3a,EA6)]; (iv) truncation of E protein by the introduction of a stop codon 9 aa 
upstream of the PBM [(A3a,E-APBM)]; (v) interruption of the E protein PBM by amino 
acid substitutions within its core sequence [(A3a,E-PBM_ )]; and (vi) replacement of the 
E protein PBM with an alternative synthetic PBM sequence with proven binding to PDZ 
domains (46) [(A3a,E-PBM*)]. The viability of each mutant was determined after rescue 
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FIG 8 Analysis of the subcellular localization of SARS-CoV 3a and E proteins by immunofluorescence. 
Vero E6 cells were grown on coverslips and infected with rSARS-CoV at a MOI of 0.3 and were 
subsequently fixed with 4% paraformaldehyde at 24 hpi. (A) Cells were labeled with antibodies specific 
for 3a protein (green) or E protein (red). (B) Cells were labeled with antibodies specific for 3a protein 
(shown in green) or for PDI (ER marker), 58 K (Golgi marker), or aconitase 2 (mitochondrion marker) 
(shown in red). In all cases, nuclei were stained with DAPI (blue). 








from infectious bacterial artificial chromosome (BAC) clones. For each mutant, two 
independently generated BACs were analyzed (Fig. 10A). Mutations within the 
N-terminal domain or TMD of protein E did not affect virus viability, indicating that 
neither harbored the motif responsible for viroporin replacement. However, when the 
protein E PBM was disrupted by amino acid substitutions [(A3a,E-PBM_)], infectious 
virus was not rescued. Consistent with this, when the original PBM was substituted for 
an alternative one [(A3a,E-PBM*)], viable virus was recovered. Furthermore, the trun- 
cated mutant [(A3a,E-APBM)] rapidly reverted in the first passage, losing the engineered 
stop codon during passage in cell culture and regaining the PBM [A3a,E-APBM(rev)]. 
Collectively, these data strongly support the idea of the requirement of the PBM for 
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FIG9 Analysis of the subcellular localization of SARS-CoV 3a by immunofluorescence. Vero E6 cells were 
grown on coverslips and infected with rSARS-CoV at a MOI of 0.3 and were fixed with 4% paraformaldehyde 
at 24 hpi. Cells were labeled with antibodies specific for 3a protein (shown in green) or for Na+/K+ ATPase 


(plasma membrane marker), rab5 (early endosome marker), rab7 (late endosome marker), and lamp-1 
(lysosome marker) (shown in red). In all cases, nuclei were stained with DAPI (blue). 


viability of SARS-CoV (Fig. 10A) and indicate that it was the E protein with a functional 
PBM that compensated for the loss of the full-length 3a protein. An additional mutant 
was generated in which the protein 3a PBM was mutated in the context of full-length 
E protein deletion [(3a-PBM~,AE)], but the rescue of this virus was not possible 
(Fig. 10B), again showing the necessity of at least one of the two proteins with a 
functional PBM for virus viability in the context of complete deletion of either the E 
gene or the 3a gene. 

Also, virus titers were significantly decreased after one passage for mutants A3a,EA1 
and A3a,E-APBM(rev) compared to the A3a virus (Fig. 10A). In addition, A3a viruses with 
small deletions throughout the C terminus of E protein were not viable, possibly due 
to the requirement of the 3a protein native structure or of a length essential for PBM 
availability for binding to other viral or cellular proteins (Fig. 10A). To analyze the 
relevance of protein shortening or folding, an additional mutant was constructed in the 
A3a background by filling in the EA4 deletion with an alanine-rich sequence in order to 
restore the original length or folding of the E protein [(A3a,EA4*)] (Fig. S2). This mutant 
virus was viable, indicating the lack of sequence specificity for the observed function- 
ality. An alternative role for other domains located in the middle of the E protein C 
terminus cannot be completely ruled out. 
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FIG 10 Mapping of the E protein domain required for the replacement of the 3a protein. (A) Viability of a set of SARS-CoV A3a mutants 
with added mutations within the E gene was evaluated in Vero E6 cells. Changes to the E protein included the following: alanine 
substitutions in the N-terminal domain [rSARS-CoV-MA15-(A3a,EA1)]; disruption of the TMD, abrogating E protein IC activity [(A3a,E- 
N15A)]; short in-frame deletions throughout the C-terminal domain [(A3a,EA2), (A3a,EA3), (A3a,EA4), (A3a,EA5), and (A3a,EA6)]; truncation 
of the E protein by introduction of a stop codon 6 aa upstream of the PBM [(A3a,E-APBM)]; interruption of the E protein PBM by amino 
acid substitutions within its core sequence [(A3a,E-PBM_)]; and replacement of the E protein PBM with an alternative PBM sequence 
[(A3a,E-PBM*)]. The APBM(rev) sequence represents the sequence of a spontaneous revertant virus isolated after growth of the original 
APBM strain in cell culture. (B) Viability of SARS-CoV AE mutants with or without 3a protein PBM was evaluated in Vero E6 cells. Titers of 
the viable viruses after the first passage (representative of three replicate experiments) are shown. Data represent means + SD. NR, not 


rescued. 





To further study replication dependence on the PBMs of E and 3a, a collection of 
mutants with one [(3a-PBM7) or (E-PBM_)] or both [(3a,E]-PBM~ ) viral PBM sequences 
removed was generated (Fig. 10A). All three PBM mutants were viable, showing that the 
consequences of the presence of E and 3a protein PBMs when these proteins have 
been completely deleted are different from those seen when only their PBM is missing. 
However, virus titers decreased 10-fold when both viral PBMs were missing in com- 
parison to the levels seen with the parental virus or mutants lacking only one PBM 
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FIG 11 Virulence of recombinant SARS-CoV combining the knockdown of 3a and E protein PBMs. (A) Schematic 
of recombinant mutants with knockdown of the PBMs of the 3a and E proteins in different combinations. Titers of 
the viable viruses (representative of three replicate experiments) are shown. Data represent means + SD. (B) 
Groups of five 16-week-old BALB/c mice were subjected to mock infection (PBS, gray filled squares) or infected with 
100,000 PFU of the parental (wild-type) virus (WT; black filled hexagons) or of genetically engineered viruses lacking 
either the protein 3a PBM (3a-PBM, red filled diamonds) or the E protein PBM (E-PBM~, ochre filled circles) or both 
[(3a,E)-PBM~, green filled circles]. Mean levels of weight loss (left) and survival (right) through 10 dpi are 
represented for each group. Data summarize two replicate experiments with equivalent results. Error bars represent 
the standard deviations of mouse weight data. 


(Fig. 11A). Remarkably, when PBM was present in E protein, the virus was highly virulent 
independently of whether the 3a protein included a PBM (Fig. 11B). In contrast, 
mortality significantly decreased when E protein lacked a PBM, regardless of the 
presence or absence of a PBM within protein 3a. These results indicated that the impact 
of the presence or absence of the E protein PBM on virulence is definitive, as it 
determines whether a virus is pathogenic or nonpathogenic, respectively. In contrast, 
the presence or absence of 3a protein with or without its PBM had little impact on virus 
virulence, illustrating the much greater relevance of the E protein PBM than of the 3a 
PBM with respect to virus pathogenicity. 


DISCUSSION 

Viroporins are highly relevant for viral replication and pathogenesis, and their 
requirement in a large number of physiological processes makes their study a field of 
growing interest (17, 18, 23). CoVs usually encode two or more viroporins, including the 
conserved structural E protein and additional ones encoded by accessory genes. The 
roles of these viroporins in replication and virulence had been studied in detail only for 
SARS-CoV E protein (34, 46, 58, 60), though roles in replication have also been 
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established for HCoV-OC43 ns12.9 (25), HCoV-229E 4a (24), and PEDV 3 (26). MERS-CoV 
encodes two proteins, E and 5, homologues of the SARS-CoV E and 3a proteins, 
respectively. While the IC activity of MERS-CoV E protein has been previously described 
(15), the potential IC activity of protein 5 is yet to be studied in detail. However, due to 
its similarity to SARS-CoV 3a, HCoV-229E 4a, and PEDV 3, IC conductance by this protein 
is expected. Furthermore, both the MERS-CoV E and 5 proteins have a putative PBM at 
their carboxy terminus, similarly to the SARS-CoV E and 3a proteins. 

SARS-CoV encodes three viroporins: 3a, E, and 8a (14, 27, 28). We have previously 
shown that a SARS-CoV mutant lacking E protein was attenuated in mice (61). Here, we 
show that removal of the 3a protein slightly reduced virus titers in vitro compared to 
the results seen with the parental virus (Fig. 1A). In contrast, the titers of the A3a mutant 
were reduced 10-fold in vivo (Fig. 1B), indicating that the 3a protein was required for 
optimal SARS-CoV replication. This result is in agreement with previous studies showing 
a slight reduction of SARS-CoV-A3a titers in Vero E6 cells (62). Nevertheless, it has to be 
noted that the mice in the previous study were infected with a human SARS-CoV Urbani 
strain, which causes only a mild murine infection, in contrast to our mouse-adapted 
strain (61). The reduction of titers after infection with A3a mutant viruses may reflect 
the role of 3a in membrane rearrangement, increasing the levels of intracellular vesicles 
that can promote nonlytic release of viral particles (30). Alternatively, protein 3a may 
induce apoptosis (30, 63) and may enhance inflammation by activating nuclear factor 
kappa B (NF-kB), leading to the production of proinflammatory cytokines such as 
interleukin-8 (IL-8) and RANTES (CCL5) (29). In fact, histopathological analysis of lungs 
from SARS-CoV-A3a-infected mice showed minimal damage or cellular infiltration at 2 
and 4 dpi, whereas mice infected with the SARS-CoV parental virus revealed interstitial 
and peribronchial cell infiltration and edema in both alveolar and bronchiolar airways 
at 2 dpi and, mainly, at 4 dpi (see Fig. S3 in the supplemental material). This suggests 
that the attenuation of the A3a mutant may be due to its inability to activate an 
exacerbated proinflammatory response, resulting in survival of infected mice. In fact, 
similar results were obtained in the lungs of mice infected with rSARS-CoV-MA15-AE, 
which also led to attenuation by downregulation of the host proinflammatory response 
(61, 64). 

We showed that deletion of SARS-CoV protein 8a alone did not have a measurable 
effect on replication and virulence in mice (Fig. 1 and 2). This is in line with the fact that 
ORF8 was lost during the SARS-CoV pandemic and that viruses lacking this gene were 
recovered from patients who died from SARS-CoV infection, supporting the idea that it 
is also not required for virulence in humans (65). However, viruses lacking 8a protein 
when 3a protein was also deleted had titers 10-fold lower than those seen with the 
parental virus, indicating that simultaneous deletion of the two proteins may contribute 
to replication to at least a small extent. 

Residues involved in E protein IC activity have been previously identified (33), and 
this activity was essential for replication and virulence (34). We previously reported that 
the E ion channel activates the inflammasome through calcium release from intracel- 
lular stores (34, 58). IC activity of SARS-CoV E protein is exerted in the membrane of the 
ERGIC. This facilitates the release of Ca?* from this intracellular compartment, which 
contributes to the activation of the inflammasome complex, leading to the release of 
proinflammatory cytokines such as tumor necrosis factor alpha (TNF-a), IL-1, and IL-6. 
The accumulation of these cytokines promotes an exacerbated proinflammatory re- 
sponse, which leads to death. Thus, E protein IC activity is a virulence factor similar to 
the IC activities of M2 from influenza virus (19) or rotavirus NSP4 (13, 66). As E protein 
IC activity promotes virulence via inflammasome activation (34, 58), we hypothesized 
that protein 3a IC activity could act in a similar way. However, as no correlation was 
found between protein 3a IC activity and SARS-CoV titer and pathogenicity, we 
conclude that, in our mouse model, the IC activity of 3a protein did not affect 
replication and virulence in the same way as that of E protein. However, SARS-CoV 3a 
protein TMD3~ and TMD[2,3]~ mutants displayed no IC activity and were attenuated. 
In principle, a possible explanation for this observation is that the TMD3~ three-amino- 
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acid mutations present in both mutant TMD3~ and mutant TMD[2,3]~ could have 
disrupted a function of the 3a protein other than its IC activity. Nevertheless, viruses 
that included the point mutations that disrupt IC activity (Y91A, H93A, and Y109A) and 
which are less likely to affect other functions of the 3a protein than the 3-amino-acid 
and 5-amino-acid mutations of TMD3~ and TMD[2,3]~, respectively, were completely 
virulent. The TMD3~ mutant showed 10-fold-lower titers in the lungs of infected mice. 
The impact of these mutants on virus replication and virulence may in principle be due 
to an exclusive effect on 3a protein functions different from its IC activity. The 
TMD[2,3]~— mutant maintained the same titers as the native virus, which could be 
explained if the mutations in TMD2 were structurally compensated for the ones in 
TMD3 when the two were simultaneously present in protein 3a. Overall, these results 
indicate that protein 3a IC activity is not involved in virus virulence. 

There are several factors that may account for the differences in the degrees of 
relevance of the SARS-CoV E and protein 3a IC activities. The two proteins localize to 
different subcellular compartments (Fig. 8A), suggesting that the IC activity of the two 
proteins regulates ion transport between different compartments. E protein locates in 
the ERGIC (56) and induces calcium efflux during SARS-CoV infection, which activates 
the inflammasome complex, leading to the acute proinflammatory response associated 
with virus pathogenicity (34, 58). As protein 3a is not found at any of the main 
intracellular Ca2* storage locations (Fig. 8B), its IC activity most likely induces cellular 
pathways that are less relevant to pathogenesis. 

Note that SARS-CoV 3a protein subcellular localization had been previously studied 
by transfecting cells with a plasmid, which overexpresses a tagged variant of 3a protein, 
leading to the conclusion that protein 3a is located at the Golgi compartment or at 
the cell surface (31, 67). These studies have also been performed in the context of the 
infection but without the use of subcellular compartment markers, leading to the 
observation that protein 3a was located in the cell membrane, in the cytoplasm, and in 
the nucleus (27). After our analysis in the context of the infection using markers for 
cellular compartments, we could not conclude that the SARS-CoV 3a protein is located 
in any of the previously described compartments. Furthermore, knowledge of the 
cellular location in which the 3a protein accumulates remained elusive after our analysis 
performed using markers for the following different intracellular compartments: ER, 
Golgi apparatus, mitochondria, early and late endosomes, lysosomes, and the plasma 
membrane. 

The E protein PBM is another virulence factor of SARS-CoV (46), with a homologous 
motif present in protein 3a. However, analysis of mutants lacking the 3a PBM showed 
no effect on virus production or virulence (Fig. 7). The requirement for protein E with 
its PBM was dominant for SARS-CoV virulence in comparison to the requirement for 
protein 3a with its PBM (Fig. 11). Furthermore, a SARS-CoV with E protein that included 
its PBM was virulent in the presence or absence of 3a protein PBM, whereas in the 
reverse situation (the presence of 3a protein with its PBM and of E protein lacking its 
PBM), the virus was always attenuated, reinforcing the idea of the dominance of the E 
protein PBM for virus pathogenicity over that of 3a protein. 

The interaction of the E protein PBM with syntenin PDZ motifs activates the p38 
mitogen-activated protein kinase (MAPK) pathway and promotes an acute proinflam- 
matory response that leads to death (46). However, other viral PBMs likely show a 
preference for distinct PDZ domains. There are 266 motifs present in more than 400 
cellular protein isoforms, each containing between 1 and 13 PDZ domains (42). As 
SARS-CoV E and 3a protein PBM sequences are different (DLLV and SVPL, respectively), 
they are likely to interact with different networks of PDZ-containing cellular proteins. 
The protein 3a PBM interaction with cellular PDZ proteins most likely induces a 
signaling pathway that either is not pathogenic for the host or is activated at a lower 
intensity due to a reduced affinity of the viral PBM for the cellular PDZ domain, 
explaining why viruses with or without the PBM of 3a showed no significant differences 
in pathogenicity. However, given that the protein 3a PBM is required for virus viability 
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when full-length E protein is missing, it can be confirmed that the pathways activated 
by the protein 3a PBM have some impact on virus replication. 

In order to study the interdependence of the SARS-CoV viroporins, the effects of 
deletion of one viroporin or of simultaneous deletions of two or three viroporins were 
determined. The 8a protein was included to determine its potential relevance in the 
absence of the other viroporins, showing a significant impact on viral growth under 
conditions in which 3a protein or E protein was absent. In contrast, we observed that 
SARS-CoV was not viable when both the 3a and E proteins were absent but that its 
viability was rescued by the presence of either 3a or E. We showed that the E protein 
could compensate for the loss of the other viroporin providing that it carried a 
functional PBM (Fig. 10A), a conclusion reinforced by the results obtained with five 
different recombinant viruses. Even more biologically relevant was the rapid reversion 
of a stop codon introduced 9 aa upstream of the carboxy terminus of the E protein in 
mutant APBM. Reversion was identified only in this mutant, presumably because 
reversion of a single altered codon occurs more easily than the 4-aa changes in E-PBM— 
or the deletions performed (Fig. 10A). Our results also showed that a 3a protein carrying 
a PBM could compensate for the loss of E protein and could restore virus replication but 
not virulence (Fig. 10B and 11). 

The idea of the requirement of SARS-CoV PBMs for virus replication and virulence is 
also supported by our previous observations revealing that AE mutants evolved to 
introduce a new transmembrane protein with a PBM to compensate for the loss of the 
whole E protein during passage either in cell culture or in vivo (68). In addition, when 
SARS-CoV-AE was passaged in mice, it spontaneously gained an internal PBM in the 8a 
protein. PBMs are also phylogenetically conserved in proteins from SARS-CoV and 
MERS-CoV isolates from humans, civet cats, dromedary camels, and bats (Fig. S4), 
further reinforcing the idea of the relevance of PBMs for coronavirus viability. Viral 
PBM-cellular PDZ interactions have been previously described in several viruses. For 
instance, the PBM of E6 oncoprotein from human papillomavirus 16 (HPV-16) interacts 
with several cellular proteins containing PDZ motifs, leading to tumorigenesis and virus 
dissemination (69). On the other hand, different strains of rabies virus (RABV) include a 
variety of PBMs in the G protein that differentially contribute to viral virulence through 
activation of a variety of signaling pathways (70). In the case of influenza virus, NS1 
protein has a PBM located at its carboxy terminus that is possibly involved in enhanced 
virulence (71). In addition, the vaccinia virus F11 protein has both a PDZ domain and 
a PBM, which influence viral spread (72). 

Overall, we conclude that the SARS-CoV E and 3a proteins have in common two 
important characteristics: IC activity and PBM. Both the IC and the PBM of E protein are 
involved in SARS-CoV virulence and replication, whereas the corresponding motifs in 
protein 3a are not. We showed a dominance of E protein PBM and IC activities over that 
of the 3a protein homologues. However, the protein 3a PBM became relevant for virus 
viability in the absence of full-length E protein. These results contribute to a better 
understanding of the role of CoV IC activities and PBMs, with impacts on the rational 
design of future vaccines and antivirals. 


MATERIALS AND METHODS 


Ethics statement. Animal experimental protocols were approved by the Ethical Committee of the 
Center for Animal Health Research (CISA-INIA) (permit numbers 2011-009 and 2011-09) in strict accor- 
dance with Spanish National Royal Decree (RD 53/2013) and international EU guidelines 2010/2063/UE 
and Spanish national law 32/2007. Infected mice were housed in a self-contained ventilated rack 
(Allentown, NJ). 

Viruses. Mouse-adapted (MA15) (48) parental wild-type (wt) and recombinant viruses were rescued 
from infectious cDNA clones generated in a bacterial artificial chromosome (BAC) (49, 73-75). 

Generation of recombinant viruses. Viruses with mutations in SARS-CoV viroporins E, 3a, and 8a 
were constructed in an infectious cDNA clone of SARS-CoV-MA15 within a bacterial artificial chromosome 
(BAC) (plasmid pBAC-SARS-CoV-MA15) (49, 73, 76). Generation of the E deletion mutant was described 
previously (61). The 3a gene was deleted by overlap extension PCR using pBAC-SARS-CoV-MA15 and the 
primers shown in Table S2 in the supplemental material. Mutations included deletion of a region 
(nt 25270 to 25668) of the SARS-CoV genome, resulting in deletion of the 3a protein while retaining the 
3b protein; a disruption of the ATG start codon of the 3a gene; and point mutations at nt 25673 and 
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25683 to introduce two stop codons and a point mutation at nt 26042 to disrupt a potential initiation 
codon. A PCR product was generated from nt 24937 to 26060 of the SARS-CoV genome, digested at 
flanking Swal and BamHI sites, and cloned into intermediate plasmid pBAC-SARS-Pmel-BamHI-SARS-CoV 
(which contains the nt-18404-to-26044 sequence of the SARS-CoV infectious cDNA clone) to generate 
plasmid pBAC-SARS-Pmel-BamHI-A3aa. Finally, that plasmid was digested with Pmel and BamHI and the 
fragment carrying the 3a deletion was reinserted into a similarly digested pBAC-SARS-CoV-MA15 plasmid 
to generate pBAC-SARS-CoV-A3a. To construct the 8a deletion mutant, a DNA fragment containing 
nt 26790 to 28753 of the SARS-CoV genome flanked by restriction sites Xcml and Nhel was assembled 
by overlap extension PCR using pBAC-SARS-CoV-MA15 as the template and the primers indicated in 
Table S1; this process resulted in the deletion of the first 82 nt of the 8a gene without affecting the 8b 
gene. The final PCR products were digested with Xcml and Nhel and cloned into intermediate plasmid 
pBAC-BamHI-RsrlI-SARS-CoV (which contains nt 26044 to 29782 of the SARS-CoV infectious cDNA clone 
[49]), thus obtaining pBAC-BamHI-Rsrll-SARS-CoV-A8a. That plasmid was digested with BamHI and Rsrll, 
and the fragment carrying the deletion of the 8a gene was reinserted into pBAC-SARS-CoV-MA15, 
generating pBAC-SARS-CoV-A8a. 

To generate mutants with mutations in the 3a IC, DNA fragments containing nt 25016 to 26044 of 
the SARS-CoV genome were produced by overlap extension PCR using plasmid pBAC-SARS-CoV-MA15 as 
a template and the primers indicated in Table S1. The following mutations were introduced: S40A (TCA 
to GCA), S48A (AGC to GCC), Y91A (TAT to GCT), H93A (CAT to GCT), Y109A (TAT to GCT), Y113A (TAT to 
GCT), Q116A (CAA to GCA), TMD1~ (including both S40A and S48A), TMD2~ (including both Y91A and 
H93A), TMD3~ (including Y109A, Y113A, and Q116A), and TMD[2,3]~ (including Y91A, H93A, Y109A, 
Y113A, and Q116A). The PCR products were digested at flanking Swal and BamHI sites and cloned into 
pBAC-Pmel-BamHI-SARS-CoV (49). Intermediate plasmids were digested with Pmel and BamHI and 
recloned into pBAC-SARS-CoV-MA15 to generate infectious SARS-CoV cDNA clones for each mutation. To 
mutate the PBM sequence of 3a, overlap extension PCR was performed using pBAC-SARS-CoV-MA15 as 
the template and the primers indicated in Table S2, resulting in a DNA fragment containing nt 26044 to 
26790 of the SARS-CoV genome. The core 3a PBM (SVPL; nt AGCGTGCCTTTG), was replaced with an 
alternative sequence (GMSM; nt GGCATGTCTATG). As the 3a and 3b genes overlap in this region, care 
was taken to introduce only silent mutations into the 3b ORF. The resulting PCR fragment was digested 
at flanking BamHI and Xcml sites and cloned into intermediate plasmid pBAC-BamHI-Rsrll-SARS-CoV (49), 
resulting in pBAC-BamHI-Rsrll-SARS-CoV-3amutPBM. This plasmid was digested with BamHI and Rsrll, and 
the fragment carrying the mutated 3a PBM was reinserted into pBAC-SARS-CoV-MA15 to generate 
plasmid pBAC-SARS-CoV-3a-PBM ~. 

Viruses combining deletion of the 3a gene with mutations throughout the E gene were generated 
by digesting previously described plasmids pBAC-SARS-CoV-MA15-EA1, -EA2, -EA3, -EA4, -EA5, and -EA6 
(75) and -E-PBM~ and -E-PBM* (46) with BamHI and Rsrll and cloning the fragments with mutations into 
a similarly digested pBAC-SARS-CoV-MA15-A3a plasmid. 

In order to generate recombinant baculoviruses for 3a protein (rBV-3a), all constructs of the SARS-CoV 
3a gene were cloned into a pFastBac vector (Invitrogen) containing a tobacco etch virus (TEV) cleavable 
site and 10 histidine residues fused to the C terminus of the 3a constructs. Recombinant baculoviruses 
were produced following the instructions of the manufacturer (Invitrogen). 

Recovery of recombinant SARS-CoV variants from the cDNA clones. BHK cells were grown to 95% 
confluence in 12.5-cm? flasks and transfected with 6 ug of infectious cDNA clone and 18 pl of 
Lipofectamine 2000 (Invitrogen), according to the manufacturer's specifications. At 6 h posttransfection 
(hpt), cells were trypsinized, added to confluent Vero E6 cells monolayers grown in 12.5-cm? flasks, and 
incubated at 37°C for 72 h. Cell supernatants were harvested and passaged once on fresh cells, and the 
recovered viruses were cloned by three rounds of plaque purification following standard procedures. 

Cells. Vero E6 cells and BHK cells were kindly provided by E. Snijder (University of Leiden, the 
Netherlands) and H. Laude (Unité de Virologie et Immunologie Molecularies, INRA, France), respectively. 
In all cases, cells were grown in Dulbecco’s modified Eagle’s medium (DMEM; Gibco) supplemented with 
25 mM HEPES, 2 mM -glutamine (Sigma), 1% nonessential amino acids (Sigma), and 10% fetal bovine 
serum (FBS; BioWhittaker, Inc.). Virus titrations were performed in Vero E6 cells as previously described 
(61). 

Mice. Eight-week-old specific-pathogen-free BALB/c Ola Hsd female mice were purchased from 
Harlan Laboratories and maintained for 8 additional weeks in the animal care facility at the National 
Center of Biotechnology (Madrid). For infection experiments, mice were anesthetized with isoflurane and 
intranasally inoculated at 16 weeks of age with 100,000 PFU of the indicated viruses. All work with 
infected animals was performed in a biosafety level 3+ (BSL3+) laboratory (CISA, INIA) by technicians 
wearing personal protection equipment (3M). 

Generation of polyclonal antibodies specific for the SARS-CoV 3a protein. A synthetic peptide 
corresponding to residues 11 to 24 of the SARS-CoV 3a protein (C-ESITAQPVKIDNAS) was generated and 
used to immunize two rabbits (Biogenes, Berlin, Germany) according to the standard protocol of the 
supplier. Serums were collected at 45 dpi and evaluated by enzyme-linked immunosorbent assay (ELISA) 
and immunofluorescence and Western blot analysis using Vero E6 cells infected with SARS-CoV-wt or 
SARS-CoV-A3a as a negative control. 

Virus genome sequencing. Regions of the SARS-CoV genome corresponding to the 3a, E, and 8a 
genes were sequenced after reverse transcriptase PCR (RT-PCR). Briefly, total RNA from infected cells or 
homogenized mouse lungs was collected and purified using an RNeasy kit (Qiagen) according to the 
manufacturer's specifications. For RT reactions, 100 ng of RNA, random oligonucleotide primers, and 
ThermoScript reverse transcriptase (Invitrogen) were used. RT products were subsequently subjected to 
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PCR using Vent polymerase (New England Biolabs) and the following primer pairs: 24937-VS (GGCGAC 
ATTTCAGGCATTAACGC) and 26086-RS (GGCACGCTAGTAGTCGTCGTCGGC), which amplify the 3a gene; 
E-VS (CTCTTCAGGAGTTGCTAATCCAGCAATGG) and E-RS (TCCAGGAGTTGTTTAAGCTTCTCAACGGTA), 
which amplify nucleotides 26017 to 26447, including the E gene; and 27545-VS (GGAGGTTCAACAAGA 
GCTCTACTCGCC) and 28008-RS (GACAGTTGATAGTAACATTAGGTGTGC), amplifying a region that includes 
the 8a gene. Sequence assembly and comparison with the parent consensus sequence were performed 
with SeqMan software (Lasergene, Madison, WI). 

Growth kinetics. Subconfluent monolayers (90% confluence) of Vero E6 cells in 12.5-cm? flasks were 
infected at a multiplicity of infection (MOI) of 0.001 with the indicated viruses. Culture supernatants were 
collected at 0, 4, 24, 48, and 72 h postinfection (hpi), and virus titers were determined as previously 
described (61). For the analysis of cell-associated virus, Vero E6 cells were infected at a MOI of 0.001 with 
the indicated viruses. At 24 and 48 hpi, cells were recovered in phosphate-buffered saline (PBS) buffer 
and disrupted by four freeze-thaw cycles. Samples were then centrifuged to remove cell debris, and 
supernatants were titrated as previously described (59). 

Virus infection and growth in mice. BALB/c mice were anesthetized with isoflurane and intranasally 
inoculated with 100,000 PFU of virus mixed with 50 wl of DMEM. Weight loss and mortality were 
evaluated daily. To determine SARS-CoV titers, lungs were homogenized in PBS containing 100 IU/ml 
penicillin, 0.1 mg/ml streptomycin, 50 ug/ml gentamicin, and 0.5 ug/ml amphotericin B (Fungizone), 
using a gentleMACS dissociator (Miltenyi Biotec, Inc.). Virus titrations were performed in Vero E6 cells as 
previously described (61). Viral titers were expressed as PFU counts per gram of tissue. 

Histopathology. Mice were sacrificed at 2 and 4 dpi. Lungs were removed, fixed in zinc formalin, and 
embedded in paraffin. Histopathological examinations were performed on sections stained with 
hematoxylin-eosin. 

lon channel reconstitution and ionic current recording. Planar bilayers were formed by apposition 
of two monolayers prepared from a mixture of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2- 
dioleoyl-sn-glycero-3-phospho-t-serine (DOPS), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine 
(DOPE) at a DOPC/DOPS/DOPE ratio of 3:1:1 (wt/wt) (Avanti Polar Lipids, Alabaster, AL) mixed in pentane 
at 5 mg/ml. Lipids were added on ~100-m-diameter orifices in the 15-ym-thick Teflon partition that 
separated two identical chambers (50, 77); the orifices were pretreated with a 1% solution of hexade- 
cane-pentane. Aqueous solutions of KCI, NaCl, or CaCl, were buffered with 5 mM HEPES at pH 6. All 
measurements were performed at room temperature (23 + 1°C). lon channel insertion was achieved by 
adding 0.5 to 1 pl of a 300 ug/ml solution of recombinant protein in a buffer containing acetonitrile- 
isopropanol (40:60) on one side of the chamber (here referred to as the cis side). 

An electric potential was applied using Ag/AgCl electrodes in 2 M KCI-1.5% agarose bridges 
assembled within standard 250-yl pipette tips. The potential was defined as positive when it was higher 
on the cis side, whereas the trans side was set to ground. An Axopatch 200B amplifier (Molecular Devices, 
Sunnyvale, CA) was used in the voltage-clamp mode to measure the current and the applied potential. 
The chamber and the head stage were isolated from external noise sources with a double metal screen 
(Amuneal Manufacturing Corp., Philadelphia, PA). Single-channel conductance data were obtained from 
current measurements under conditions of an applied potential of +100 mV and were evaluated using 
the Gaussian fit tool of Sigma plot 12 (Systat Software, Inc.). 

The reversal potential, Eev was determined as follows. First, a lipid membrane was formed at a given 
salt concentration gradient. Second, one or several channels were inserted into the bilayer and a net ionic 
current appeared due to the concentration gradient. Third, the ionic current through the channel or 
channels was manually set to zero by adjusting the applied potential. The potential needed to achieve 
zero current was then corrected using values corresponding to the liquid junction potentials of the 
electrode salt bridges (78) to obtain the final Eev: 

Confocal microscopy. Vero E6 cells were grown to 90% confluence on glass coverslips and infected 
with wt SARS-CoV at a MOI of 0.3 PFU/cell. At 24 hpi, the medium was removed and cells were washed 
twice with PBS and fixed with 4% paraformaldehyde-PBS for 30 min at room temperature. Cells were 
subsequently washed twice with PBS, permeabilized for 10 min with ice-cold methanol, and then blocked 
with PBS containing 10% FBS for 40 min at room temperature. Immunofluorescence was performed 
using mouse Abs specific for SARS-CoV E protein (generated as described in reference 56) (1:500), protein 
disulfide isomerase (PDI; Abcam, Inc.) (1:500), 58 K (Abcam, Inc.) (1:100), aconitase 2 (Abcam, Inc.) (1:500), 
rab5 (BD Biosciences) (1:100), rab7 (BD Biosciences) (1:100), lamp-1 (Santa Cruz Biotechnologies) (1:50), 
and Nat/K* ATPase (Santa Cruz Biotechnologies) (1:50) and rabbit Abs specific for SARS-CoV 3a protein 
(1:500). Primary antibodies were diluted in PBS containing 5% FBS and incubated for 90 min at room 
temperature, and then coverslips were washed four times with PBS before incubation with secondary 
antibodies was performed. Alexa 488- or Alexa 546-conjugated antibodies specific for the different 
species (Invitrogen) were diluted 1:500 in PBS containing 5% FBS and incubated for 45 min at room 
temperature. Nuclei were stained using DAPI (4’,6-diamidino-2-phenylindole; Sigma) (1:200), and cov- 
erslips were mounted in ProLong Gold anti-fade reagent (Invitrogen) and examined on a Leica SP5 
confocal microscope (Leica Microsystems, Inc.). Image analysis was performed using ImageJ (79) and the 
JACoP plug-in (57). Three areas per image of 120 by 120 pixels with high accumulation of 3a protein were 
analyzed to determine Pearson’s coefficient (Pc). Pc values below 0.6 were considered negative for 
colocalization. Partial colocalization was considered to have occurred for Pc values between 0.6 and 0.85. 
Pc values between 0.85 and 1 were considered positive for colocalization. 

Production and purification of SARS-CoV 3a protein and its mutant variants from a recombi- 
nant baculovirus (rBV-3a). H5 cells at 80% confluence were infected (MOI = 1) with rBV-3a (constructed 
as specified in the “Generation of recombinant viruses” section above) and incubated at 22°C for 72 h. 
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Cells were harvested and resuspended in lysis buffer 1 (50 mM Tris-HCl, 300 mM NaCl, 0.5% Triton X-100, 
pH 7.5) supplemented with 1% protease inhibitor cocktail (Sigma). Protein extracts were centrifuged at 
12,000 X g for 10 min at 4°C, and the pellets were resuspended in lysis buffer 2 (8 M urea, 50 mM Tris-HCl, 
300 mM NaCl, 1% IGEPAL, 1 mM B-mercaptoethanol, 10 mM imidazole, pH 7.5). Each sample was 
sonicated three times for 20 s each time and centrifuged at 12,000 X g for 10 min at 4°C. Finally, the 
protein present in the supernatant was purified through metal affinity chromatography (IMAC) using 
cobalt resin (Clontech) following the manufacturer's instructions. Every fraction from the purification 
process was analyzed in 12% polyacrylamide gels using Coomassie Blue EZBlue gel staining reagent 
(Sigma). Then, the obtained protein was desalted using a PD-10 desalting column (GE Healthcare) and 
eluted in PBS. 

Statistical analysis. Two-tailed, unpaired Student's t tests were used to analyze the differences in 
mean values between groups. All results were expressed as means + standard deviations; P values of 
<0.1 were considered significant. 


SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at https://doi.org/10.1128/mBio 
.02325-17. 
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ACKNOWLEDGMENTS 

We thank Marga Gonzalez (CNB-CSIC) for her technical assistance. In vivo experi- 
ments were performed at INIA-CISA (Madrid, Spain). 

This work was supported by grants from the Government of Spain (BIO2013-42869-R 
and BlO2016-75549-R AEI/FEDER, UE), the European Zoonotic Anticipation and Pre- 
paredness Initiative (ZAPI) (IMI_JU_115760), and the U.S. National Institutes of Health 
(NIH) (0258-3413/HHSN266200700010C awarded to L.E., 2P01AI060699 awarded to L.E. 
and S.P., and RO1 Al129269 awarded to S.P.). V.M.A. and M.Q.M. are grateful for the 
support of the Government of Spain (FIS2013-40473-P and FIS2016-75257-P AEI/FEDER, 
UE) and Universitat Jaume | (P1.1B2015-28). C.C.R. received a contract from Fundacion 
La Caixa. The funders had no role in study design, data collection and analysis, decision 
to publish, or preparation of the manuscript. 


REFERENCES 
1. Perlman S, Netland J. 2009. Coronaviruses post-SARS: update on repli- 


mBio’ 


of circulating bat coronaviruses shows potential for human emergence. 


cation and pathogenesis. Nat Rev Microbiol 7:439- 450. https://doi.org/ 
10.1038/nrmicro2147. 

. Rota PA, Oberste MS, Monroe SS, Nix WA, Campagnoli R, Icenogle JP, 
Peñaranda S, Bankamp B, Maher K, Chen MH, Tong S, Tamin A, Lowe 
L, Frace M, DeRisi JL, Chen Q, Wang D, Erdman DD, Peret TCT, Burns 
C, Ksiazek TG, Rollin PE, Sanchez A, Liffick S, Holloway B, Limor J, 
McCaustland K, Olsen-Rasmussen M, Fouchier R, Günther S, Oster- 
haus ADME, Drosten C, Pallansch MA, Anderson LJ, Bellini WJ. 2003. 
Characterization of a novel coronavirus associated with severe acute 
respiratory syndrome. Science 300:1394-1399. https://doi.org/10 
.1126/science.1085952. 

. Drosten C, Günther S, Preiser W, van der Werf S, Brodt HR, Becker S, 
Rabenau H, Panning M, Kolesnikova L, Fouchier RA, Berger A, Burguière 
AM, Cinatl J, Eickmann M, Escriou N, Grywna K, Kramme S, Manuguerra 
JC, Müller S, Rickerts V, Stürmer M, Vieth S, Klenk HD, Osterhaus AD, 
Schmitz H, Doerr HW. 2003. Identification of a novel coronavirus in 
patients with severe acute respiratory syndrome. N Engl J Med 348: 
1967-1976. https://doi.org/10.1056/NEJMoa030747. 

. Lau SK, Woo PC, Li KS, Huang Y, Tsoi HW, Wong BH, Wong SS, Leung 
SY, Chan KH, Yuen KY. 2005. Severe acute respiratory syndrome 
coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci U 
S A 102:14040-14045. https://doi.org/10.1073/pnas.0506735102. 

. Menachery VD, Yount BL, Jr, Debbink K, Agnihothram S, Gralinski LE, 
Plante JA, Graham RL, Scobey T, Ge XY, Donaldson EF, Randell SH, 
Lanzavecchia A, Marasco WA, Shi ZL, Baric RS. 2015. A SARS-like cluster 


May/June 2018 Volume 9 Issue 3 e02325-17 


Nat Med 21:1508-1513. https://doi.org/10.1038/nm.3985. 


. Menachery VD, Yount BL, Jr, Sims AC, Debbink K, Agnihothram SS, 


Gralinski LE, Graham RL, Scobey T, Plante JA, Royal SR, Swanstrom J, 
Sheahan TP, Pickles RJ, Corti D, Randell SH, Lanzavecchia A, Marasco 
WA, Baric RS. 2016. SARS-like WIV1-CoV poised for human emergence. 
Proc Natl Acad Sci U S A 113:3048-3053. https://doi.org/10.1073/pnas 
1517719113; 


. Zaki AM, van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. 


2012. Isolation of a novel coronavirus from a man with pneumonia in 
Saudi Arabia. N Engl J Med 367:1814-1820. https://doi.org/10.1056/ 
NEJMoa1211721. 


. Bermingham A, Chand M, Brown C, Aarons E, Tong C, Langrish C, 


Hoschler K, Brown K, Galiano M, Myers R, Pebody R, Green H, Boddington 
N, Gopal R, Price N, Newsholme W, Drosten C, Fouchier R, Zambon M. 
2012. Severe respiratory illness caused by a novel coronavirus, in a patient 
transferred to the United Kingdom from the Middle East, September 2012. 
Euro Surveill 17:20290. 


. Nieto-Torres JL, Verdiá-Báguena C, Castaño-Rodriguez C, Aguilella VM, 


Enjuanes L. 2015. Relevance of viroporin ion channel activity on viral 
replication and pathogenesis. Viruses 7:3552-3573. https://doi.org/10 
.3390/v7072786. 


. Pinto LH, Holsinger LJ, Lamb RA. 1992. Influenza virus M2 protein has ion 


channel activity. Cell 69:517-528. https://doi.org/10.1016/0092-8674 
(92)90452-1. 


. Ewart GD, Sutherland T, Gage PW, Cox GB. 1996. The Vpu protein of 


mbio.asm.org 21 


Hio'wseoiqu Aq pausiqnd - 810Z ‘pz Ae uo Hio"wseolqu wo papeojuMoq 


Castano-Rodriguez et al. 


12; 


13. 


14. 


15: 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25, 


26. 


27. 


28. 


29; 


30. 


May/June 2018 Volume 9 


human immunodeficiency virus type 1 forms cation-selective ion chan- 
nels. J Virol 70:7108-7115. 

Pavlović D, Neville DC, Argaud O, Blumberg B, Dwek RA, Fischer WB, 
Zitzmann N. 2003. The hepatitis C virus p7 protein forms an ion 
channel that is inhibited by long-alkyl-chain iminosugar derivatives. 
Proc Natl Acad Sci U S A 100:6104-6108. https://doi.org/10.1073/ 
pnas.1031527100. 

Hyser JM, Collinson-Pautz MR, Utama B, Estes MK. 2010. Rotavirus dis- 
rupts calcium homeostasis by NSP4 viroporin activity. MBio 1:e00265-10. 
https://doi.org/10.1128/mBio.00265-10. 

Liao Y, Lescar J, Tam JP, Liu DX. 2004. Expression of SARS-coronavirus 
envelope protein in Escherichia coli cells alters membrane permeability. 
Biochem Biophys Res Commun 325:374-380. https://doi.org/10.1016/j 
-bbrc.2004.10.050. 

Surya W, Li Y, Verdia-Baguena C, Aguilella VM, Torres J. 2015. MERS 
coronavirus envelope protein has a single transmembrane domain that 
forms pentameric ion channels. Virus Res 201:61-66. https://doi.org/10 
.1016/j.virusres.2015.02.023. 

Ao D, Sun SQ, Guo HC. 2014. Topology and biological function of 
enterovirus non-structural protein 2B as a member of the viroporin 
family. Vet Res 45:87. https://doi.org/10.1186/s13567-014-0087-6. 
Nieva JL, Madan V, Carrasco L. 2012. Viroporins: structure and biologi- 
cal functions. Nat Rev Microbiol 10:563-574. https://doi.org/10.1038/ 
nrmicro2820. 

Hover S, Foster B, Barr JN, Mankouri J. 2017. Viral dependence on cellular 
ion channels—an emerging anti-viral target? J Gen Virol 98:345-351. 
https://doi.org/10.1099/jgv.0.000712. 

Ichinohe T, Pang IK, Iwasaki A. 2010. Influenza virus activates inflam- 
masomes via its intracellular M2 ion channel. Nat Immunol 11:404-410. 
https://doi.org/10.1038/ni.1861. 

Zhang K, Hou Q, Zhong Z, Li X, Chen H, Li W, Wen J, Wang L, Liu W, 
Zhong F. 2013. Porcine reproductive and respiratory syndrome virus 
activates inflammasomes of porcine alveolar macrophages via its small 
envelope protein E. Virology 442:156-162. https://doi.org/10.1016/j.virol 
.2013.04.007. 

Ito M, Yanagi Y, Ichinohe T. 2012. Encephalomyocarditis virus viroporin 
2B activates NLRP3 inflammasome. PLoS Pathog 8:e1002857. https://doi 
.org/10.1371/journal.ppat.1002857. 

Andrew A, Strebel K. 2010. HIV-1 Vpu targets cell surface markers CD4 
and BST-2 through distinct mechanisms. Mol Aspects Med 31:407-417. 
https://doi.org/10.1016/j.mam.2010.08.002. 

Griffin SD. 2009. Plugging the holes in hepatitis C virus antiviral therapy. 
Proc Natl Acad Sci U S A 106:12567-12568. https://doi.org/10.1073/pnas 
.0906760106. 

Zhang R, Wang K, Lv W, Yu W, Xie S, Xu K, Schwarz W, Xiong S, Sun B. 
2014. The ORF4a protein of human coronavirus 229E functions as a 
viroporin that regulates viral production. Biochim Biophys Acta 1838: 
1088-1095. https://doi.org/10.1016/j.bbamem.2013.07.025. 

Zhang R, Wang K, Ping X, Yu W, Qian Z, Xiong S, Sun B. 2015. The ns12.9 
accessory protein of human coronavirus OC43 is a viroporin involved in 
virion morphogenesis and pathogenesis. J Virol 89:11383-11395. 
https://doi.org/10.1128/JVI.01986-15. 

Wang K, Lu W, Chen J, Xie S, Shi H, Hsu H, Yu W, Xu K, Bian C, Fischer WB, 
Schwarz W, Feng L, Sun B. 2012. PEDV ORF3 encodes an ion channel 
protein and regulates virus production. FEBS Lett 586:384-391. https:// 
doi.org/10.1016/j.febslet.2012.01.005. 

Lu W, Zheng BJ, Xu K, Schwarz W, Du L, Wong CK, Chen J, Duan S, 
Deubel V, Sun B. 2006. Severe acute respiratory syndrome-associated 
coronavirus 3a protein forms an ion channel and modulates virus re- 
lease. Proc Natl Acad Sci U S A 103:12540-12545. https://doi.org/10 
.1073/pnas.0605402103. 

Chen CC, Krüger J, Sramala |, Hsu HJ, Henklein P, Chen YM, Fischer WB. 
2011. ORF8a of SARS-CoV forms an ion channel: experiments and mo- 
lecular dynamics simulations. Biochim Biophys Acta 1808:572-579. 
https://doi.org/10.1016/j.bbamem.2010.08.004. 

Kanzawa N, Nishigaki K, Hayashi T, Ishii Y, Furukawa S, Niiro A, Yasui F, 
Kohara M, Morita K, Matsushima K, Le MQ, Masuda T, Kannagi M. 2006. 
Augmentation of chemokine production by severe acute respiratory 
syndrome coronavirus 3a/X1 and 7a/X4 proteins through NF-kappaB 
activation. FEBS Lett 580:6807-6812. https://doi.org/10.1016/j.febslet 
.2006.11.046. 

Freundt EC, Yu L, Goldsmith CS, Welsh S, Cheng A, Yount B, Liu W, 
Frieman MB, Buchholz UJ, Screaton GR, Lippincott-Schwartz J, Zaki SR, 
Xu XN, Baric RS, Subbarao K, Lenardo MJ. 2010. The open reading frame 


Issue 3 e02325-17 


3k 


32. 


33% 


34. 


35, 


36. 


37, 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


48. 


49. 


mBio’ 


3a protein of severe acute respiratory syndrome-associated coronavirus 
promotes membrane rearrangement and cell death. J Virol 84: 
1097-1109. https://doi.org/10.1128/JVI.01662-09. 

Yuan X, Li J, Shan Y, Yang Z, Zhao Z, Chen B, Yao Z, Dong B, Wang S, 
Chen J, Cong Y. 2005. Subcellular localization and membrane association 
of SARS-CoV 3a protein. Virus Res 109:191-202. https://doi.org/10.1016/ 
j.virusres.2005.01.001. 

Vennema H, Godeke GJ, Rossen JWA, Voorhout WF, Horzinek MC, Op- 
stelten DJ, Rottier PJM. 1996. Nucleocapsid-independent assembly of 
coronavirus-like particles by co-expression of viral envelope protein 
genes. EMBO J 15:2020-2028. 

Verdia-Baguena C, Nieto-Torres JL, Alcaraz A, Dediego ML, Torres J, 
Aguilella VM, Enjuanes L. 2012. Coronavirus E protein forms ion channels 
with functionally and structurally involved membrane lipids. Virology 
432:485- 494. https://doi.org/10.1016/j.virol.2012.07.005. 

Nieto-Torres JL, Dediego ML, Verdia-Baguena C, Jimenez-Guardeno JM, 
Regla-Nava JA, Fernandez-Delgado R, Castafo-Rodriguez C, Alcaraz A, 
Torres J, Aguilella VM, Enjuanes L. 2014. Severe acute respiratory syn- 
drome coronavirus envelope protein ion channel activity promotes virus 
fitness and pathogenesis. PLoS Pathog 10:e1004077. https://doi.org/10 
.1371/journal.ppat.1004077. 

Oostra M, de Haan CA, Rottier PJ. 2007. The 29-nucleotide deletion 
present in human but not in animal severe acute respiratory syndrome 
coronaviruses disrupts the functional expression of open reading frame 
8. J Virol 81:13876 -13888. https://doi.org/10.1128/JVI.01631-07. 

Chen CY, Ping YH, Lee HC, Chen KH, Lee YM, Chan YJ, Lien TC, Jap TS, 
Lin CH, Kao LS, Chen YM. 2007. Open reading frame 8a of the human 
severe acute respiratory syndrome coronavirus not only promotes viral 
replication but also induces apoptosis. J Infect Dis 196:405- 415. https:// 
doi.org/10.1086/519166. 

Piserchio A, Spaller M, Mierke DF. 2006. Targeting the PDZ domains of 
molecular scaffolds of transmembrane ion channels. AAPS J 
8:E396 -E401. https://doi.org/10.1007/BF02854911. 

Feng W, Zhang M. 2009. Organization and dynamics of PDZ-domain- 
related supramodules in the postsynaptic density. Nat Rev Neurosci 
10:87-99. https://doi.org/10.1038/nrn2540. 

Guggino WB, Stanton BA. 2006. New insights into cystic fibrosis: molec- 
ular switches that regulate CFTR. Nat Rev Mol Cell Biol 7:426-436. 
https://doi.org/10.1038/nrm1949. 

Hung AY, Sheng M. 2002. PDZ domains: structural modules for protein 
complex assembly. J Biol Chem 277:5699 -5702. https://doi.org/10.1074/ 
jbc.R100065200. 

Münz M, Hein J, Biggin PC. 2012. The role of flexibility and conforma- 
tional selection in the binding promiscuity of PDZ domains. PLoS Com- 
put Biol 8:e1002749. https://doi.org/10.1371/journal.pcbi.1002749. 
Luck K, Charbonnier S, Travé G. 2012. The emerging contribution of 
sequence context to the specificity of protein interactions mediated 
by PDZ domains. FEBS Lett 586:2648-2661. https://doi.org/10.1016/j 
febslet.2012.03.056. 

Javier RT, Rice AP. 2011. Emerging theme: cellular PDZ proteins as 
common targets of pathogenic viruses. J Virol 85:11544-11556. https:// 
doi.org/10.1128/JVI.05410-11. 

Lenfant N, Polanowska J, Bamps S, Omi S, Borg JP, Reboul J. 2010. A 
genome-wide study of PDZ-domain interactions in C. elegans reveals a 
high frequency of non-canonical binding. BMC Genomics 11:671. 
https://doi.org/10.1186/1471-2164-11-671. 

Gallardo R, Ivarsson Y, Schymkowitz J, Rousseau F, Zimmermann P. 2010. 
Structural diversity of PDZ-lipid interactions. ChemBioChem 11: 
456-467. https://doi.org/10.1002/cbic.200900616. 

Jimenez-Guardefto JM, Nieto-Torres JL, DeDiego ML, Regla-Nava JA, 
Fernandez-Delgado R, Castaho-Rodriguez C, Enjuanes L. 2014. The PDZ- 
binding motif of severe acute respiratory syndrome coronavirus enve- 
lope protein is a determinant of viral pathogenesis. PLoS Pathog 10: 
e1004320. https://doi.org/10.1371/journal.ppat.1004320. 

Teoh KT, Siu YL, Chan WL, Schlüter MA, Liu CJ, Peiris JS, Bruzzone R, 
Margolis B, Nal B. 2010. The SARS coronavirus E protein interacts with 
PALS1 and alters tight junction formation and epithelial morphogenesis. 
Mol Biol Cell 21:3838-3852. https://doi.org/10.1091/mbc.E10-04-0338. 
Roberts A, Deming D, Paddock CD, Cheng A, Yount B, Vogel L, Herman 
BD, Sheahan T, Heise M, Genrich GL, Zaki SR, Baric R, Subbarao K. 2007. 
A mouse-adapted SARS-coronavirus causes disease and mortality in 
BALB/c mice. PLoS Pathog 3:e5. https://doi.org/10.1371/journal.ppat 
0030005. 

Almazan F, DeDiego ML, Galan C, Escors D, Alvarez E, Ortego J, Sola |, 


mbio.asm.org 22 


Hio'wseoiqu Aq pausiqnd - 810Z ‘pz Ae uo Hio"wseolqu wo papeojuMoq 


SARS-CoV Viroporins in Replication and Pathogenesis 


50. 


51. 


52. 


53; 


54. 


59; 


56. 


57. 


58. 


59, 


60. 


61. 


62. 


63. 


64. 


65. 


May/June 2018 Volume 9 


Zuñiga S, Alonso S, Moreno JL, Nogales A, Capiscol C, Enjuanes L. 2006. 
Construction of a severe acute respiratory syndrome coronavirus infec- 
tious cDNA clone and a replicon to study coronavirus RNA synthesis. J 
Virol 80:10900-10906. https://doi.org/10.1128/JVI.00385-06. 

Bezrukov SM, Vodyanoy |. 1993. Probing alamethicin channels with 
water-soluble polymers. Effect on conductance of channel states. Bio- 
phys J 64:16-25. https://doi.org/10.1016/S0006-3495(93)81336-5. 
Aguilella VM, Queralt-Martín M, Aguilella-Arzo M, Alcaraz A. 2011. In- 
sights on the permeability of wide protein channels: measurement and 
interpretation of ion selectivity. Integr Biol 3:159-172. https://doi.org/ 
10.1039/c0ib00048e. 

Hille B. 2001. lon channels of excitable membranes, 3rd ed. Sinauer 
Associates, Sunderland, MA. 

Hodgkin AL, Katz B. 1949. The effect of sodium ions on the electrical 
activity of giant axon of the squid. J Physiol 108:37-77. https://doi.org/ 
10.1113/jphysiol.1949.sp004310. 

Hsu HJ, Fischer WB. 2012. In silico investigations of possible routes of 
assembly of ORF 3a from SARS-CoV. J Mol Model 18:501-514. https:// 
doi.org/10.1007/s00894-011-1092-6. 

Krüger J, Fischer WB. 2009. Assembly of viral membrane proteins. J 
Chem Theory Comput 5:2503-2513. https://doi.org/10.1021/ct900185n. 
Nieto-Torres JL, Dediego ML, Alvarez E, Jiménez-Guardefo JM, Regla- 
Nava JA, Llorente M, Kremer L, Shuo S, Enjuanes L. 2011. Subcellular 
location and topology of severe acute respiratory syndrome coronavirus 
envelope protein. Virology 415:69-82. https://doi.org/10.1016/j.virol 
.2011.03.029. 

Bolte S, Cordeliéres FP. 2006. A guided tour into subcellular colocaliza- 
tion analysis in light microscopy. J Microsc 224:213-232. https://doi.org/ 
10.1111/j.1365-2818.2006.01706.x. 

Nieto-Torres JL, Verdia-Baguena C, Jimenez-Guardeho JM, Regla-Nava 
JA, Castafio-Rodriguez C, Fernandez-Delgado R, Torres J, Aguilella VM, 
Enjuanes L. 2015. Severe acute respiratory syndrome coronavirus E 
protein transports calcium ions and activates the NLRP3 inflammasome. 
Virology 485:330-339. https://doi.org/10.1016/j.virol.2015.08.010. 
Berridge MJ, Lipp P, Bootman MD. 2000. The versatility and universality 
of calcium signalling. Nat Rev Mol Cell Biol 1:11-21. https://doi.org/10 
.1038/35036035. 

DeDiego ML, Nieto-Torres JL, Jimenez-Guardefo JM, Regla-Nava JA, 
Alvarez E, Oliveros JC, Zhao J, Fett C, Perlman S, Enjuanes L. 2011. Severe 
acute respiratory syndrome coronavirus envelope protein regulates cell 
stress response and apoptosis. PLoS Pathog 7:e1002315. https://doi.org/ 
10.1371/journal.ppat.1002315. 

DeDiego ML, Alvarez E, Almazan F, Rejas MT, Lamirande E, Roberts A, 
Shieh WJ, Zaki SR, Subbarao K, Enjuanes L. 2007. A severe acute respi- 
ratory syndrome coronavirus that lacks the E gene is attenuated in vitro 
and in vivo. J Virol 81:1701-1713. https://doi.org/10.1128/JVI.01467-06. 
Yount B, Roberts RS, Sims AC, Deming D, Frieman MB, Sparks J, Denison 
MR, Davis N, Baric RS. 2005. Severe acute respiratory syndrome corona- 
virus group-specific open reading frames encode nonessential functions 
for replication in cell cultures and mice. J Virol 79:14909-14922. https:// 
doi.org/10.1128/JVI.79.23.14909-14922.2005. 

Chan CM, Tsoi H, Chan WM, Zhai S, Wong CO, Yao X, Chan WY, Tsui SK, 
Chan HY. 2009. The ion channel activity of the SARS-coronavirus 3a 
protein is linked to its pro-apoptotic function. Int J Biochem Cell Biol 
41:2232-2239. https://doi.org/10.1016/j.biocel.2009.04.019. 

DeDiego ML, Nieto-Torres JL, Regla-Nava JA, Jimenez-Guardeno JM, 
Fernandez-Delgado R, Fett C, Castaho-Rodriguez C, Perlman S, Enjuanes 
L. 2014. Inhibition of NF-kB-mediated inflammation in severe acute 
respiratory syndrome coronavirus-infected mice increases survival. J 
Virol 88:913-924. https://doi.org/10.1128/JVI.02576-13. 

Chiu RW, Chim SS, Tong YK, Fung KS, Chan PK, Zhao GP, Lo YM. 2005. 


Issue 3 e02325-17 


66. 


67. 


68. 


69. 


70. 


71. 


72. 


73. 


74. 


75. 


76. 


77. 


78. 


79. 


mBio’ 


Tracing SARS-coronavirus variant with large genomic deletion. Emerg 
Infect Dis 11:168-170. https://doi.org/10.3201/eid1101.040544. 

Zhang M, Zeng CQ, Morris AP, Estes MK. 2000. A functional NSP4 
enterotoxin peptide secreted from rotavirus-infected cells. J Virol 74: 
11663-11670. https://doi.org/10.1128/JVI.74.24.11663-11670.2000. 
Minakshi R, Padhan K. 2014. The YXXPhi motif within the severe acute 
respiratory syndrome coronavirus (SARS-CoV) 3a protein is crucial for its 
intracellular transport. Virol J 11:75. https://doi.org/10.1186/1743-422X 
211-75. 

Jimenez-Guardeho JM, Regla-Nava JA, Nieto-Torres JL, DeDiego ML, 
Castaho-Rodriguez C, Fernandez-Delgado R, Perlman S, Enjuanes L. 2015. 
Identification of the mechanisms causing reversion to virulence in an 
attenuated SARS-CoV for the design of a genetically stable vaccine. PLoS 
Pathog 11:e1005215. https://doi.org/10.1371/journal.ppat.1005215. 

Ganti K, Broniarczyk J, Manoubi W, Massimi P, Mittal S, Pim D, Szalmas A, 
Thatte J, Thomas M, Tomaić V, Banks L. 2015. The human papillomavirus 
E6 PDZ binding motif: from life cycle to malignancy. Viruses 
7:3530-3551. https://doi.org/10.3390/v7072785. 

Caillet-Saguy C, Maisonneuve P, Delhommel F, Terrien E, Babault N, 
Lafon M, Cordier F, Wolff N. 2015. Strategies to interfere with PDZ- 
mediated interactions in neurons: what we can learn from the rabies 
virus. Prog Biophys Mol Biol 119:53-59. https://doi.org/10.1016/j 
.pbiomolbio.2015.02.007. 

Jackson D, Hossain MJ, Hickman D, Perez DR, Lamb RA. 2008. A new 
influenza virus virulence determinant: the NS1 protein four C-terminal 
residues modulate pathogenicity. Proc Natl Acad Sci U S A 105: 
4381-4386. https://doi.org/10.1073/pnas.0800482105. 

Handa Y, Durkin CH, Dodding MP, Way M. 2013. Vaccinia virus F11 
promotes viral spread by acting as a PDZ-containing scaffolding protein 
to bind myosin-9A and inhibit RhoA signaling. Cell Host Microbe 14: 
51-62. https://doi.org/10.1016/j.chom.2013.06.006. 

Fett C, DeDiego ML, Regla-Nava JA, Enjuanes L, Perlman S. 2013. Com- 
plete protection against severe acute respiratory syndrome coronavirus- 
mediated lethal respiratory disease in aged mice by immunization with 
a mouse-adapted virus lacking E protein. J Virol 87:6551-6559. https:// 
doi.org/10.1128/JVI.00087-13. 

DeDiego ML, Pewe L, Alvarez E, Rejas MT, Perlman S, Enjuanes L. 2008. 
Pathogenicity of severe acute respiratory coronavirus deletion mutants 
in hACE-2 transgenic mice. Virology 376:379-389. https://doi.org/10 
.1016/j.virol.2008.03.005. 

Regla-Nava JA, Nieto-Torres JL, Jimenez-Guardeho JM, Fernandez- 
Delgado R, Fett C, Castafo-Rodriguez C, Perlman S, Enjuanes L, DeDiego 
ML. 2015. Severe acute respiratory syndrome coronaviruses with muta- 
tions in the E protein are attenuated and promising vaccine candidates. 
J Virol 89:3870 -3887. https://doi.org/10.1128/JVI.03566-14. 

DeDiego ML, Nieto-Torres JL, Jimenez-Guardeño JM, Regla-Nava JA, 
Castaño-Rodriguez C, Fernandez-Delgado R, Usera F, Enjuanes L. 2014. 
Coronavirus virulence genes with main focus on SARS-CoV envelope 
gene. Virus Res 194:124-137. https://doi.org/10.1016/j.virusres.2014.07 
.024. 

Montal M, Mueller P. 1972. Formation of bimolecular membranes from 
lipid monolayers and a study of their electrical properties. Proc Natl Acad 
Sci U S A 69:3561-3566. https://doi.org/10.1073/pnas.69.12.3561. 
Alcaraz A, Nestorovich EM, López ML, García-Giménez E, Bezrukov SM, 
Aguilella VM. 2009. Diffusion, exclusion, and specific binding in a large 
channel: a study of OmpF selectivity inversion. Biophys J 96:56-66. 
https://doi.org/10.1016/j.bpj.2008.09.024. 

Schneider CA, Rasband WS, Eliceiri KW. 2012. NIH Image to ImageJ: 25 
years of image analysis. Nat Methods 9:671-675. https://doi.org/10 
.1038/nmeth.2089. 


mbio.asm.org 23 


Hio' wise oiqu Aq pausiqnd - 810Z ‘pz Aen uo Bio'wse'oqu wo papeojuMoq 


